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ABSTRACT 


A  frequency  dependent  analysis  of  a  shielded  edge-guided 
mode  isolator  is  presented.  A  Fourier  transform  technique  is 
applied  to  the  boundary  expressions  of  a  structure  built  on 
a  dielectric  substrate,  and  the  resulting  equations  are 
solved  for  the  wavelength  ratio.  By  using  perturbation 
analysis  and  the  results  obtained  for  the  dielectric  case, 
solutions  for  the  normalized  propagation  constant  and 
attenuation  for  waves  traveling  in  the  -Z  and  +Z  directions, 
in  a  structure  built  on  a  ferrite  substrate,  are  obtained. 
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I .  INTRODUCTION 


The  analysis  of  various  configurations  of  microstrip 
transmission  lines,  is  of  great  importance,  primarily,  to 
the  industrial  area,  since  these  kinds  of  transmission  lines 
are  easy  to  manufacture,  and  they  are  suitable  for  use  in 
microwave  integrated  circuits  due  to  their  small  dimensions. 

Two  major  groups  of  microstrip  transmission  lines  can 
be  defined,  related  to  the  substrate  material  that  is  used. 
These  are  the  dielectric  and  ferrite  substrate  groups.  This 
study  presents  the  analysis  of  a  new  type  of  ferrite  built, 
shielded  microstrip  isolator,  using  the  edge-guided  mode  of 
propaga ti on . 

It  was  shown  by  M.  E.  Hines  [Ref.  1]  that  an  edge-guided 
mode  of  propagation  occurs  in  a  wide  microstrip  transmission 
line  using  a  ferrite  slab  magnetized  perpendicular  to  the 
ground  plane.  The  R.  F.  fields  patterns  in  such  a  structure 
are  shown  in  Figure  1  [Ref.  1]. 

Propagation  occurs  in  both  j^Z  directions  with  equal 
phase  velocities  and  loss,  and  the  R.  F.  fields  patterns  are 
mirror  images  for  both  directions  of  propagation.  In  the 
dominant  mode  the  energy  is  shifted  from  one  side  to  the 
other  with  reversal  of  the  direction  of  propagation.  Hines 
suggested  that  this  phenomenon  could  be  used  in  the  analysis 
of  nonreciprocal  devices  such  as  isolators,  phase  shifters 
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X 


Figure  1.  R.F.  Fields  of  Dominant  Mode  in  Micro- 
strip  Line.  D.C.  Magnetic  Field  in 
-y  Direction. 
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and  other  components,  where  nonreciprocal  behavior  can  be 
attained  by  perturbing  the  structure  in  an  asymmetrical  way. 

To  achieve  isolation,  a  resistive  load  can  be  placed 
only  on  one  side  of  the  structure;  thus  high  transmission 
loss  occurs  for  the  propagation  direction  for  which  the 
energy  is  concentrated  on  the  lossy  side,  and  lower  loss 
results  for  the  opposite  direction  where  the  energy  is 
concentrated  on  the  other  side.  The  disadvantage  of  this 
field  displacement  isolator  was  the  small  ratio  between 
backward  and  forward  losses,  which  was  not  sufficient  for 
isolation  purposes. 

An  improved  edge-guided  mode  isolator  was  suggested  by 
K.  Araki ,  T.  Koyama  and  Y.  Naito  [Ref.  2]  for  which  no  lossy 
electric  material  was  used,  but  instead  one  edge  of  the 
conducting  strip  was  shorted  to  ground. 

This  isolation  structure,  shown  in  Figure  2,  was  built 
and  tested  and  the  experimental  result  showed  large  attenua¬ 
tion  in  the  backward  direction  and  small  insertion  loss  in 
the  forward  direction. 

The  following  study  is  an  investigation  of  the  theo¬ 
retical  behavior  of  a  model  for  this  isolator  structure. 

The  exact  configuration  of  the  isolator  model  is  shown 
in  Figure  3.  Since  most  of  the  field  is  confined  in  the 
dielectric  substrate  then  the  left,  right  and  upper  conducting 
walls  of  the  shield  have  negligible  affect  on  the  electric 
and  magnetic  fields  configuration,  thus  the  model  which  is 
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Figure  2.  Isolator  Structure. 

easier  to  analyze  can  replace  the  open  boundary-isolator 
configuration  [Figure  2].  The  shield  is  a  rectangular  wave¬ 
guide  which  implies  that  the  highest  frequency  of  operation 
should  not  exceed  the  cut-off  frequency  of  the  mode  of 

the  waveguide. 

The  mathematical  method  which  was  selected  for  this 
study  was  used  in  earlier  works  like  references  3,  4,  and 
5.  This  mathematical  method  can  be  used  for  calculating 
wavelength  ratio  and  characteristic  impedance  for  non  TEM, 

TE  or  TM  transmission  lines.  Essentially,  this  method  deals 
with  the  boundary  conditions  of  a  structure,  after  they  were 
transformed  into  the  spectral  domain.  Then  by  using  the 


10 


method  of  moment  [Ref.  4]  and  assuming  either  current  density 
distributions  or  electric  field  distributions  (depending  on 
the  structure)  numerical  solutions  for  the  wavelength  ratio 
and  the  characteristic  impedance  can  be  obtained. 

In  essence,  the  analysis  procedure  is  a  complete 
solution  for  the  dispersion  characteristics  of  a  structure 
built  on  a  dielectric  substrate,  and  then  by  applying 
perturbation  theory,  the  phase  and  the  attenuation  constants 
of  a  structure  built  on  a  ferrite  substrate  are  obtained. 

This  ferrite  built  device  is  the  isolator. 


REGION  1  (AIR) 

^  \i/  ■ 

H 

^  W  «-  1 

r  egion  '  2 

T 

0 

Figure  3.  Three  dimensional  Isolator  and  Regions. 
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II .  DISPERSION  CHARACTERISTICS  ON  DIELECTRIC  SUBSTRATE 


A.  FIELD  AND  BOUNDARY  CONDITIONS 

Assume  waves  are  propagating  along  the  structure  shown 
in  Figure  3,  in  the  -Z  direction  such  that  the  propagation 
term  is  where  m  is  the  radian  frequency,  and  Y  =  ci+j3 

is  the  propagation  constant.  We  assume  from  now  on  that  for 
the  dielectric  substrate  y  is  pure  imaginary,  hence  a=0. 

Later  we  shall  see  that  for  the  ferrite  subststrate,  y  is  a 
complex  number  including  both  terms  a  and  3. 

The  longitudinal  components  of  the  electric  and  magnetic 

fields  can  be  expressed  in  terms  of  the  scalar  electric 

e  h 

potential  and  the  scalar  magnetic  potential  cj)  as 


II 

rsi 

UJ 

Kc.2 

(x,y)e^^ 

(1) 

N 

II 

kc.2 

(x,y)e^^ 

(2) 

where  Kc^^  =  K^. 

i  =  1.2 

defining  the 

spatial 

regions  of  the 

Structure  (Figure 

3).  Using  the 

two 

f i  el  d 

expressions  (1) 

and  (2) 

together 

wi th  Maxwe 1 1 ' s 

curl 

equations  leads 

to  the 

transverse 

components  of 

the 

electric 

and  magnetic  fields  as 

f ol 1 ows : 
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Ex^ 

3^- 

=  (y—  - 

34' 

(3) 

Hx. 

3^; 

3y  ^ 

(4) 

3^? 

=  (y^  +  juv. 

3  (j)^ 

T>r^ 

(5) 

Hy,. 

34,^ 

=  (y—  - 

3<J>® 

(6) 

i  =  1 

,2. 

Applyi ng 

boundary  conditions 

at  the 

interface  between 

the  dielectric 

substrate  and 

the 

bottom 

wall  of  the 

waveguide 

(y=0),  and  assuming  perfect 

conducting  walls,  tangential 

electric  fields  should  be  zero  or 

expl i ci t ly  wri ti ng 

Bx2 

(x.o.z)  =  0 

(7) 

EZ2 

(x,0,z)  =  0. 

(8) 

At  the  interface  between  regions  1 

and  2  ,  (y=D) 

the 

tangential  electric  field  must  be 

con ti nuous  : 

Ex^ 

(x.D.z)  =  EX2 

(x.D, 

z) 

(9) 

(x,D,z)  =  EZ2 

(x  ,D, 

z). 

(10) 

The  electric  fields  at  the  interface  between  regions 
1  and  2  exist  only  in  the  dielectric  portion  of  the  interface 
and  can  be  written  as: 
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Ex^ (x ,D  ,z  ) 


on  strip 
(x)e^^ 


W<X<A 


/  0  on  s  tri p 
ie2(x)e^^  W<X<A. 


(n) 


(12) 


The  tangential  magnetic  fields  at  the  interface  are 
discontinuous  by  the  surface  current  densities,  and  assuming 
i nf i ni tess i mal ly  thin  strip  (can  be  considered  as  current 
sheet)  the  following  boundary  conditions  can  be  written  at 
y  =  D 

rO  W<X<A 

Hx,  (x.D.z)  -  Hx„(x,D,z)  =  \ 

vjz(x)e^  on  strip 


(0  W<X<A  .... 

Hz, {x,D,z)  -  Hz«(x,D,z)  =  ^  ^  ^ 

(jx(x)e^  on  strip. 

At  the  interface  between  region  1  and  the  upper  wall  of 
the  waveguide  (y=H)  the  tangential  electric  fields  must  be 
zero  or  explicitly  writing 


Ex^{x,H,z)  =  0 

(15) 

Ez^  (x,H,z)  =  0. 

(16) 

Substituting  the  field  expressions  of  equations  (1) 
through  (6)  into  the  boundary  condition  expressions  of 
equations  (7)  through  (16)  yields  the  following  equations: 
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3't'2(x.0) 

3X 


jup2 


l(x.O) 


3y 


=  0 


(17) 


KC2^  '|)|(x,0)  =  0 


(18) 


a^®(x.D)  a*T(x,D) 


3*5(X,D) 


3x 


Juv 


1  3y 


- ^“>^2 - W 


■(>®(X,D)  =  KC2^<(i|(X.D) 


(19) 

(20) 


a<t®(x,D) 

r - 33J - j“yi 


3i)>!J(X,D)  (0  on  strip 


ay 


ex(x)  W<X<A 


(21) 


Kc 


T  <(>®(X.D)  = 


0  on  strip 
e^Cx)  W<X<A 


3*"(X.D)  a*®(X.D) 

'  +  jcoel  ' 


3X 


ay 


(22) 


(23) 


-  (y 


a^2(x>D) 


3x 


+  j  ue2 


9<^o(X,D) 


3y 


-)  = 


ro  w<x<A 

)jz(x)  on  strip 


^i,‘(X,D)  -  Kc^  <j)2(X,D)  = 


0  W<X<A 
Jx(x)  on  strip 


(X,H) 


8X 


-  j(^Pl 


ai>!j(x,H) 

ay 


=  0 


$®(x.H)  =  0. 


(24) 

(25) 

(26) 
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B.  SPECTRAL  DOMAIN  TRANSFORM 

0  h 

The  scalar  potential  functions  c|)^.  and  must  satisfy 
Helmholtz's  equations  in  the  two  spatial  regions,  thus 

Kc?)  <f®(x.y)  =  0  (27) 

+  Kc?)  <(>'?(x.y)  =  0  (28) 

2  2  2  2  2  2 
where  Kc^.  =  y  +  K^.  =  K^.  -  B  and  vxy  is  a  two-dimensional 

Laplacian  operator. 

As  was  suggested  by  Itoh  and  Mittra  [Ref.  3]  a  Fourier 
transform  was  used  in  the  a  domain.  The  transform  is  defined 
by: 

00 

(x,y)]  =  'J'^-(a,y)  =  /  <t.^.(x,y)  dx  (29) 

i  =  1.2. 

Since  the  structure  is  closed  boundary,  a  finite  Fourier 
transform  must  be  used  instead  of  the  infinite  Fourier 
transform.  The  finite  Fourier  transform  is  given  by 

a 

F^[*,-(x,y)]  =  <t^.(ccny)  =  /  dx  (30) 

0 

.  2  •  TT  •  n 

where  an  =  — — 

and 

9  ^  ( X  ,  y ) 

FxC— T3; - ^  Fx^^-(x,y)].  (31) 
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By  using  equations  (30)  and  (31),  a  general  transform 
of  equations  (27)  and  (28)  is  given  by 


3(ti?(x,y)  '  3'(i?(x,y)  p 

J  ^  KcfF^E-),.  (x,y)]  =  0  (32) 

9  X  ^ 


-1 

3y 


Or  explicitly  writing: 


a  <5.(ar\y)  2 

- ^ ^^(an,y) 

9y 


(33) 


where  y?  =  +  6^  -  K.^  i  =  1,2. 


Equation  (33)  should  be  analyzed  for  both  regions  1  and 
2,  and  the  solutions  will  be  the  transforms  of  the  electric 
and  magnetic  scalar  potentials  for  both  regions. 

For  region  1 , 

+  6^  -  (34) 

Where  and  B  =  ^  . 

Substituting  and  a  into  equation  (34)  obtains, 

^  (35) 
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X'  is  the  effective  structure  wavelength  and  is  related  to 
the  free  space  wavelength  X  by 


xo 


>  X ' 


XQ 


(36) 


where  is  the  effective  relative  dielectric  constant 

of  the  device. 

Since  e  r-j  =  1  in  region  1,  equation  (36)  can  be 
rewritten  as 


xo>x'  =  — -  .  (37) 

■''erTrf 

By  using  equation  (37)  together  with  equation  (35)  one  can 
find  that 

>  CL  ^  (38) 


Therefore  is  always  a  real  quantity. 
For  region  2: 


2 

^2  = 


2 


(39) 


2 

where  K2  =  w  • 

By  substituting  <2  and  S  into  equation  (39)  one  can  obtain 
’2'  •  ■  (|?)']  .  (40) 
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From  equation  (40)  it  is  clear  that  can  be  either 
a  real  or  an  imaginary  quantity,  depending  on  the  value  of 
an . 

Y2  will  be  imaginary  for 


and  wi 1 1  be  real  for 


.  .  <  Ur2^r2  ' 


2n  J  /  Xo  X “  <  an  < 

ly-  - 


AO 


^ur,.r2  - 


(42) 


The  last  two  equations  should  be  carefully  observed  while 
solving  equation  (33). 

For  region  1  equation  (33)  can  be  written  as 
2 


3$,  (an,y)  2 

- P -  =  Yi  <^»^(an,y) 

ay 

and  the  solution  has  the  form  of 

<J»l(an,y)  =  A(an)cosh  yi*y+  B(an)  sinh  y^^-y 


(43) 


(44) 


the 


For  region  2,  two  solutions  do  exist  corresponding  to 
real  and  imaginary  values  of  y2* 
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For  Y2  > 

’  C(an)  cosh  Y2y  +  D(an)  sinh 
and  for  y imaginary 

II  II 

<^’2('^'^>y)”C(an)cosY2  y  +  jD(an)  sin  Y2  y  (^^6) 

II 

where  Y2  =  JY2 

After  knowing  the  solutions  for  ^-jCanjy)  and 

e  e  h  h 

one  can  write  the  transforms  of  (j)^  ,  (^2^  *  ^2  ^  both  regions 

as  follows. 


Region  1  D£y<_H 

^e 


^’^(an,y)  =  A  (an ) coshY*!  (y-D )  +  B  (an)  sinh  Yi(y”D) 


(47) 


tl(oin.y)  =  A'’(an)  coshY,  (y-D)  +  B^’Can)  sinh  y,  (y-D)  (48) 


Region  2 


0£y_<D 

C^(an)  coshY2y  +  D®(an)  sinhY2y  .  (y2  i^eal  )  (49) 

e  »  e  " 

Cy{an)  COSY2  y  +  jDi-(oin)  sinY2  y.  (^2  iniaginary) 

(50) 

c[](on)  coshY2y  +  D|](an)  sinh  Y2y.  (^2  (5^) 

h  “  h  " 

Cj(an)  C0SY2  y  +  jDj(an)  sinY2  ^"^2  i”i2i9inary) 

(52) 
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Superscript  (e)  indicates  the  electric  field  case  and 
(h)  indicates  the  magnetic  field  case. 

All  coefficients  that  appear  in  equations  (47)  through 
(52)  can  be  determined  by  the  boundary  condition  expressions. 
The  above  can  be  done  by  taking  the  Fourier  transform  of 
equations  (17)  through  (26),  substituting  the  field 
expressions  for  both  cases  of  y2»  and  imaginary,  and 

finally  solving  for  the  coefficients. 

The  following  equations  are  obtained: 

Hyperbolic  Case  (y2  ■  real) 

-  joinYC^(cin  )  -  jujy2Y2DH ^ 


Kc2^C^(an)  =  0 


(54) 


KCi^A^(an)  =  ff^(an) 


(55) 


(56) 


KCi^A^(an)  -  Kc2  [cj^(an)  coshY2D  D[|(an)sinhY2D]  ■ 


=^j^(cin) 


(57) 
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(58) 


-janYA^(an)  +  juci  Y-j  B®  (a  n)  - 
h  h 

-  { -jany  [C[]  (otn )  coshy^D  +  dJJ  (an  )  si  nhy2D]  + 

+  ^  y  2^^  ( an  )  cos  hy2D]  }  =  «^^(an) 

Kc^A®(cin)  =  Kc2[C^(an)coshY2D  +  D^(an)sinhY2D]  (59) 

-  joinYA®(an)  -  juy  ^  y  i  ( a )  =  (60) 

-  jany[C^(an)coshy2D  +  D^(an)sinhy2D]  - 

h  h 

-  ^ ^  ^  ^ ^  0^  ( an )  cos hy 2^] 


"  jany[A^(an)coshy-j(H-D)  +  B^(an)sinhy-j(H-D)]  -  (61) 

-  ja)UlY-j[A^(an)sinhy^(H-D)  +  B^(an)coshy-j(H-D)]  =  0 

Kc^[A®(ctn)coshY^  (H-D)  +  B®  (a )  s  i  nhY^  ( H-0 )  ]  =  0.  (62) 


The  coefficients  for  the  hyperbolic  case  are  obtained 
as  f ol 1 ows  : 


C^(an)  =  0 

p  E.ian) 

A®(an)  =  - 

Kc^ 


(63) 

(64) 


D[](cin)  =  0 


(65) 


B*’(on)  =  -  — ^  S^(an)  +  gy(an) 


wu^y ^  KC^ 


wy-jy-j  X 


(66) 
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1 


(67) 


Du(ctn)  =  - 5-i -  S-(an) 

"  KC2‘^sinhY2D 


C[|{cin)  =  - 


any 


B®(an)  =  - 


Kc^  (jop2T2^ ^ 

ctghy-j  (H-D) 


^z(“n)  ■*■ 


5^(an) 

(68) 


Kc. 


^ -  -^2  ^ 


(69) 


u  an  •  Y  •  ctghy, ( H-D )  CtghYi(H-D) 

A"(an)  =  - ‘—2 -  E^(an)  -  j . ^ -  ?.(an)  (70) 


wy-jY  1  KC-j 


wy-j  Y-j 


When  substituting  equation  (63)  through  (70)  into 
equations  (57)  and  (58)  one  can  form  two  sets  of  equations 


F!j'{an,8)ff^(an)  +  F2”  ( an  .  6  ) 


F,^(an,e) E^(an)  +  F , ^ ( an , 8 ) S  ( on )  =  J  (an) 


(71) 

(72) 


or  in  matrix  form 


F^*^(an  ,6) 

(an  ,  B  ) 

(  “"  ) 

‘7^(an) 

F3*^(an.8) 

F^^ ( an  ,  B  ) 

^^(an) 

/  (an ) 

L.  Z  ^ 

(73) 


where  y  =  jB. 


The  elements  of  the  matrix  [F^]  are  given  as  follows: 
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F^^(an,6)  =  -  j 


Kc^^ctghy^ (H-D)  KC2^ctghY2D 


l»U-|  Y-| 


(74) 


F«  (an, 3)  =  j' 


anBctghy^ (H-D)  a^3ctghY2D 


(jy.  Y 


in 


0)112^2 


F2*^(an,3)  =  -  F2*^(an,3) 


(75) 

(76) 


F-*^(an,3)  =  j 


(an3)^ctghYi(H-D)  aj£^YiCtghYi(H-D) 


tjjyiYi  KC] 


Kc. 


-  (77) 


2 

(jiz 2^ 2^  (an3)  ctghY2D 


Kc, 


Kc 


2  ^v^2'Y2 


The  same  method  can  be  applied  to  the  trigonometric  case 
(y2  imaginary)  and  the  following  equations  are  obtained: 

c®(an)  =  0  (78) 

olfCan)  =  0  (79) 

DT(an)  =  -  j  - ! -  ff^(an)  (80) 

'  Kc2  sinY2"D 

c!|l(an)  =  - s — (an)  -  j  - - ^^(on) 

(81) 


Kc^^A^(an)  -  KC2^  [  cif  ( an  )  cos  y  2  "D  +  j  d!,!  ( an  )  s  i  n  Y2 "  D]  = 

=  ( an )  ( 82 ) 
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(83) 


-janYA*^(an)  +  jwe^YiB^(an)  -  j  anYCj  ( an  )  COSY2  "  D 


+  a3e2Y2  "^7  ( )  COSY2  "D  =  * 

When  substituting  equations  (69),  (70)  and  (78)  through 
(81)  into  equations  (82)  and  (83)  one  can  obtain  the  matrix 
form 


pj ( an  ,  3  )  F2 ( an  ,  3  ) 
(an  ,  3  )  F^ (an  ,  3  ) 


J 

£-^(011) 

(84) 


where  again  y  = 

and  the  elements  of  the  matrix  [f"^]  are  given  as  follows: 


F^  (an,3)  =  -  j 


F^  (an , 3 )  =  j 


Kc^^ctghy^ (H-D)  Kc2^ctgY2"D 


wy-iY-j  ^^2^2 

an3ctghY-j  ( H-D )  an3ctgY2''D 


wPlY-j 

F2^(an,3)  =  -  F2^(otn,3) 


ajp2Y2 


F^  (an,  3)  =  j 


(an3)^ctghYi(H"D)  a)e^Y-|CtghYi(H-D) 


w  y  1 Y 1  Kc  ^ ' 


Kc 


(85) 

(86) 

(87) 

(88) 


(an3  )‘^ctgY2"D  ims  2"  ctgy 

KC2^ 


26 


Expressions  (73)  and  (84)  can  be  written  in  a  general 


form  as 


F3’''’(c<n,8)  .  F^’’'’(c<n,8) 


r  /  xi 

Jy^{o.n) 

5’,  (an ) 

e7_  ( an ) 

L  z ' 

L  Z  '  'j 

(89) 


from  which  solutions  to  ^^(ctn)  and  ^^(an)  in  terms  of 
and  /^(an)  can  be  obtained  as  follows: 

^■^(an )  F4^  ’  ‘  (otn  ,3  )  -  *7^  ( an  )  F2^ ’"^  ( an  , 3  ) 
F^“’T(an,6)*F4H.T(^^^g)  .  F2” ’^  ( ctn  ,  6  )  •  F3^ ’^  ( an  ,  8 ) 

-  *7^  (  an  )  F^^ ’^  (  an  ,  3  )  +  «7^  (  an  )  •  Fij^ ’"^  (  an  ,  3  ) 

(91  ) 


(90) 


Define  the  following  terms 


DN  =  F'j'’''’(an,6)  F”’^(an,6)  -  F^ ’^  ( an  ,  8  )  F^ ’^  ( an  ,  6  )  (92) 


p  H.T 

,,H,T  _  '^4(an,8) 

'^Kan.S)  ■  DN 

(93) 

p  H.T 

.,H,T  _  ■  '^2(an,8) 

'^2(an,8)  ■  DN 

(94) 

F 

iijHjT  _  3(anj3)  _  M  p1 

'^3(an,6)  ■  DN  "2 

(95) 

T  .  F/>T(an,8) 

4(an,6)  DN 

(96) 

27 


which  lead  to  the  final  expressions 


E^(an)  =  ’"^(an  ,  S  )  J^ian)  +  M2*^(an,B)  ^^(an)  (97) 

=  M2’'*'(an,3)  J ^{an)  +  M^’^(an,e) 

C.  DETERMINENTAL  EQUATION 

Nothing  seems  to  be  gained  so  far  since  neither  the 
electric  field  and  the  current  density  components  nor  their 
transforms  are  known.  In  order  to  simplify  equation  (97)  and 
(98),  the  method  of  moments  [Ref.  4]  is  applied  in  the 
spectral  domain.  Define  a  scalar  product  over  the  domain 
-oo<an<«  according  to  reference  6  as 

n  =  co 

<a (an ) , b (an ) >  =  ^  a(an)*b  (an).  (99) 

n  =  -00 

Thus  when  applying  this  concept  to  equation  (95)  and  (96) 
and  choosing  a(an)  =  £'^(an)  or  j^(a.r\)  and  b(an)  =  e/^(an)  or 
o/^(an)  respectively,  one  can  obtain 

<M!j* ’'*'(an  ,  6  )  ^/^(an),  J^(an)>  +  <M2 (an  ,$  )^2  (an  )  ,  j^(an  )>  =  0 

(100) 

<M2’'*'(an  ,6)o7^(an)  ,J'2(an)>  +  ’*^  (  an  ,  3  )  (  an  )  ,  (  an  )  >  =  0. 

(101) 
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By  using  Parseval's  theorem  it  can  easily  be  shown  that 
the  right  hand  side  of  the  above  two  equations  is  equal  to 
zero,  due  to  the  orthogonality  of  S^{a.r\)  and  ^^(an),  and 
S^^(an)  and  *f^(an) . 

In  general,  one  can  expand  each  current  density  component 
in  a  set  of  basis  functions  such  as: 

oo 

02(x)  =  f^^(x)  -  (102) 

I<=1 

00 

k=l 

In  the  following  analysis,  one  term  approximation  is  used 
since  it  was  found  to  be  less  complicated,  without  much 
degradation  in  the  accuracy  of  the  final  results. 

So  with 

J^ian)  '  F{al«f^(x)}  =  a-jf^Cctn)  (104) 

and 

J-^(an)  =  F{a2-g^(x)}  =  a2-gj^(oin)  (105) 

equations  (100)  and  (101)  become 

00  oo 

^2  S  ’"^(an  ,  6  )  I  g^  (an)  I  ^  ^1  X]  M2’^(an,6) 

n= -«  n= 

(106) 

f^{an)g*^{an)  =  0 
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(107) 


^  Y'  Ii^^’^(ctn,6)g^{ctn)  f*(an)  +  a-j  ^M5’^(an,6)' 

n  =  -00 

•|fz(“n)|2  =  0. 

From  the  geometry  of  the  device  (Fig.  3)  and  the  loca¬ 
tion  of  the  coordinate  system  it  is  clear  that  all  choices  of 
current  density  distributions,  and  J^(x),  are  neither 

even  nor  odd  functions. 

Each  current  density  distribution  can  be  expressed  as 
a  linear  function  of  even  and  odd  functions  as  follows: 

Jz(x)  =  "l^^ze^^)  ^  ^zo 

^  (x)]  (109) 

When  taking  the  Fourier  transform  of  equations  (108) 
and  (109)  then  according  to  reference  7,  the  following 
transforms  are  obtained 


•/^(ocn)  =  a^[f^^(an) 

(no) 

+  jgj^o(cin)]  =  a2gj^(an). 

(in) 

In  order  to  solve  equations  (106)  and  (107)  for  the 
non-trivial  solution,  the  determinant  of  the  coefficient 
matrix  should  be  set  to  zero  for  all  sets  of  physical  para¬ 
meters,  at  each  frequency  of  operation.  The  above  can  be 
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achieved  by  finding  the  appropriate  3  that  satisfies  this 
requirement.  From  that  3  the  dispersion  characteristics  of 
the  device  can  be  calculated. 

The  determi nental  equation  has  the  form  of 


M!j*’'^(an  ,3)  i  g^(an)  |  ^ 

n  =  -00 


M^’^(an,3) I f^Can) 

n  =  -oo 


^  M2’^(an,3)  •  f^ian)  g*(an) 

n=  -00 


M^’'^(an,3) 

_n  =  -oo 


g)^(an)fz(an] 


=  0 


(112) 


By  investigating  equations  (93)  through  (96)  one  may 
observe  that  M-j(an,3)  and  M^(an,3)  are  even  functions,  while 
M2(cxn,3)  and  M2(ctn,3)  are  odd  functions. 

Using  this  information  together  with  equations  (110)  and 
(111)  leads  to  the  final  expression  of  the  determi nental 
equati on 


00  ■" 
(an  ,3  )  1  gx(  )  1  ^ 

• 

(an  ,3  )  1  f^  (an  )  i  ^ 

_n  =  -oo 

-n  =  -‘»  - 

M2’'^(an,3)  •  B(an) 


Ln  =  - 


=  0 


(113) 


where 

B(an)  =  -  f2g(an)  •  9x0^*^^^  (H^) 
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In  addition  to  the  requirement  that  the  determinant  of 
the  coefficients  matrix  should  be  set  to  zero,  one  can  see 
from  equations  (106),  (113)  and  (114)  that  the  ratio  of  the 


two  coefficients  has  the  form  of 


(115) 


M2’^(an,3)  B(an) 


n  =  -00 


Since  the  expression  in  the  brackets  is  a  real  quantity, 
it  is  well  understood  that  there  is  a  90°  phase  difference 
(in  time)  between  the  two  current  density  components. 

D.  CURRENT  DENSITY  COMPONENTS 

After  the  determi nental  equation  was  simplified  as  shown 
in  equation  (113),  the  two  current  density  components  v3^(x) 
and  'lj^(x)  were  approximated,  thus  that  the  equation  could 
have  been  solved.  Various  one  term  approximations  were 
investigated  and  the  set  of  components  that  was  chosen  to  be 
substituted  in  the  determi nental  equation  was  the  one  found 
as  the  best  approximation  of  both  current  distribution 
components . 

The  approximated  current  density  component  in  the 
z-direction  has  the  form 


on  the  strip 


el sewhere 


(116) 
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and  in  the  x-di recti  on 


5^2  cos^  on  the  strip 

(117) 

0  elsewhere 

The  two  current  density  components  are  shown  in  Figure  4. 


Figure  4.  Longitudinal  and  Transverse 
Current  Density  Components. 

It  is  clear  that  the  ratio  ^3/^1  =  =  4.85‘10^ 

indicates  that  most  of  the  current  that  flows  in  the  z-direction 
is  concentrated  at  the  right  edge  of  the  strip  [Fig.  3].  On 
the  other  hand  the  choice  of  indicates  that  there  is  no 

current  in  the  x-direction  at  the  right-edge  of  the  strip. 
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Computer  output  indicated  a  frequency  dependent  of  the 


ratio  a^/  as  shown  in  Figure  5.  This  curve  was  plotted  for 

D  =0.125"  w  =  0.45"  =  16  and  a,  =  10'^°. 

r  I 


FREQ  (GHz) 


Figure  5.  Ratio  of  the  Maximum  Current 
Density  Components 

From  Figure  5  one  may  learn  that  the  transverse  current 
component  on  the  strip  increases  relative  to  the  longitudinal 
component,  as  the  frequency  increases. 
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The  Fourier  transforms  of  the  current  density  components 
are  gi ven  as  : 


J  ^{ar\) 


ai*f^(an)  = 


(118) 


^(e^^cosanw-1  )  +  e^^ •  an  •  s i  nanw 

^ - 2 - 

(^)  +  (an) 


+  j 


20  2020 
-an(e  cosanw-1)  +  — e  sinanw 
. 

(|^)  +  (an) 


J’^(an) 


-  /  ,  ^janw  ^ 

-  j  "2‘=°^2^7  " 


®  =  a2.g^(an)  = 


(119) 


,  5-“  COSanw  TT  . 

=  a^  2w  sinanw  +  an 

2  - 2 - 2  +  J  2w _ 

_(k)  -  (-)  (|^)'  -  (an)' 

where  the  even  and  odd  components  of  each  transform  can  be 
identified  by  refering  to  equations  (110)  and  (111). 

After  substituting  equations  (118)  and  (119)  into 
equation  (113),  the  determi nent al  equation  was  programmed  into 
Fortran  language  and  the  root  3  that  set  the  determi nenta 1 
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equation  equal  to  zero,  was  found  by  an  iteration  method  on 
a  digital  computer.  From  the  value  of  6,  X‘/X  could  easily 
be  calculated. 

Figure  6  shows: 

1.  Computed  and  measured  X'/X  vs.  frequency  for 

D  =0.125“,  W  =0.45“,  =  16. 

2.  Computed  X'/\  vs.  frequency  for  D  =0.025", 

w  =0.45",  =  16. 

Both  curves  are  plotted  for  A  =  0.9",  H  =  0.4". 

It  can  be  seen  that  there  is  a  good  agreement  between 
experimental  and  computed  wavelength  ratio.  The  highest 
deviation  is  on  the  order  of  10  percent  which  can  be  attributed 
largely  to  the  experimental  apparatus. 

One  can  learn  that  the  wavelength's  ratio  is  directly 
proportional  to  the  dielectric  substrate  width',  which  implies 
that  the  fields  are  more  confined  in  the  dielectric  substrate 
as  its  width  decreases. 

E.  AVERAGE  POWER  FLOW 

For  later  discussion  of  the  characteristics  of  the  device 
built  on  ferrite  substrate,  a  perturbation  technique  is  used. 
One  of  the  entries  to  the  perturbation  expression  is  the 
average  power  flow  in  the  device  built  on  a  dielectric 
substrate . 

A  general  expression  for  the  time  average  power  flow  in 
the  direction  of  propagation  is  given  by 
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SHIELD  DIMENSIONS- 


0.28- 


0.27H - r 

1  2 


3 


FREQ  CGHz^ 


Fiqure  6,  Comouted  and  Measured  Dispersion  Characteristics 
for  D  =  0.025"  W  =  0.45"  and  D  =  0.125" 

W  =  0.45"  with  £  =16. 

r 
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(120) 


Pave  =  i  Re  //  (ExH  )»a_  da 
^  S  ^ 

Since  (ExH  )*a^  =  ExHy  -  EyHx* 

and  da  =  dx  dy 

then  equation  (120)  can  be  rewritten  as 

Pave  =  -  i  Re  //  (ExHy*  -  EyHx*)  dx  dy  (121) 

S 

The  minus  sign  was  added  in  equation  (121)  to  make  the 
result  positive  (recall  propagation  in  the  -z  direction  was 
assumed ) . 

Both  the  electric  and  the  magnetic  fields  are  known  in 

e  h 

terms  of  the  scalar  potentials  ^  and  (p  thus  substituting 
equations  (3)  through  (6)  into  equation  (121)  and  using 
Y  =  JB  for  the  waves  traveling  in  the  -z  direction,  leads  to 
the  following  expression 


P^ ave  =-2^6 


d  (p  i  ^  d  cp 

^ 9x  "  ay 


ay 


ff  34.^  3<Ci  a*®  I  ,  , 

^  dxdyj  (122) 

i  =  l  .2 

Subscript  i  denotes  the  spatial  regions,  air  and  dielectric 
substrate  respectively,  for  which  equation  (122)  must  be 
eval uated . 
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The  total  average  power  flow  in  the  device  is  the  sum 
of  the  spatial  power  components  that  flow  both  in  the  air  and 
in  the  dielectric  substrate  as  shown  in  the  following  equation 

PaveioT  "  ^1^^'® 

Applying  Parseval's  theorem  to  equation  (122)  obtains 


P.ave  =  -  ^  Re 


n  =  -oo  Ja 


an  oi)3e^  |  (an  ,y  )  | 


(124) 


-  an  a)0y  I  $^.  ( an  ,  3  )  I  - 


(an,y) 


9y 


-  cj  3  e 


9(j)?(an,y)  ^ 
9y  I 


jan3 


^?(an.y) - ^ 


(an,y) 


h*  9^^(an,y) 


ay 


+  j  a  n  K 


*i 


3<i’J(cin,y) 


3y 


^h* 

^  ^  ^ - , 


dy 


The  limits  of  integration  depend  on  each  region  and  are 
H  D 

for  region  1  and  J'  for  region  2. 

D  d 

After  obtaining  the  general  expression  for  average  power 
flow,  one  can  apply  it  for  both  regions. 

Recall  equations  (47),  (48),  (64),  (66),  (69)  and  (70) 
and  have  the  first  two  modified,  the  following  equations  can 
be  obtained 
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sin  hYi{y-H) 


(125) 


^^(an.y)  =  A  (an) 


'^^(an.y)  =  A  (an)  cos  hYi(y-H) 


(126) 


where 


A  (an)  = 


S’,  (an) 


Kc^'^sin  hy^  (D-H) 


A  (an)  =  -  j 


anSS^ (an ) 
Kci^ 


(an ) 


(127) 

(128) 


(DV^ Yl  sin  h  (D-H) 


After  using  the  modified  expressions  for  $®(an,y)  and 
*^(an,y)  in  region  1,  and  integrating  with  respect  to  y, 
one  can  obtain 


P^ave  =  Re 


(H-D) 


E 


'0)6 


sin  h2Y^ (H-D) 


2y 


P  -e  2  p  -h  ; 
an  e-i|A  (an)|  +  Yi  U-||A  (an)| 


+  ujS 


sinh2YT (H-D) 

27^ 


+  (H-P) 


-h 


-e 


+  Y-j  I  A  (an  )  I 


+  jan 


an  y 1  I  A  ( an  )  I  + 
e  -h* 


6  A  (an)A  (an) 


-  K^'^A  (an)A  (an) 


sinhZy^ (H-D) 


(129) 
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In  region  2,  there  are  two  expressions  depending  upon  whether 
Y2  is  a  real  or  imaginary  quantity.  As  in  region  1,  the 
integration  with  respect  to  y  can  be  done  analytically,  and 
the  following  equations  are  obtained. 

For  real  , 

00 

p^aven  =  {a  Re 

n  =  -00 

+  -^2  I  c[j(an  )  I  ^ 


’si  nh2Y2D 


2yc 


-  D 


+  uB 


2 

an  Uo 


I  Cm (an ) I ^  + 


0)  3 


2 

an  Zr 


I D^(an  )  I ' 


(130) 


+  Y2^^2  I  ^ 


■ 

"sinh2YoD 

r 

- ^  +  D 

""2 

+  jan 

3^D^(an)cj^  (an) 


K2^D®*(an)  c[](an) 


sinh  2Y2D 
and  for  Y2  imaginary,  (y2  =  jY2") 


n  =  -00 


D  - 


si  n2Y2"D’ 


2yo" 


+  0)  3 


n2  irh|  + 

an  Ii2  I  I  ^ 


II 2  [  e  I  2 

+  Y  2  ^  2  ‘  ^ 


D  + 


si n2Y2  "D 


2y," 


(131) 


+  an 


2  e  h*  2  h 
3  DjCl|!  +  K^^Cjdj 


sin2Y2  D 
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D^(ctn),  c|j(cin),  Dj(an)  and  Cj(an)  are  given  in  equation  (67), 
(68),  (80)  and  (81)  respectively 

The  total  power  in  region  2  becomes 

P2ave  ^  (132)- 

Detailed  development  of  equations  (129)  through  (131) 
is  given  in  appendix  A. 

Curves  of  a-domain  power  distributions  for  regions  1  and 
2  are  shown  in  Figures  7  and  8  respectively.  The  ratio  of  the 
power  carried  in  each  region  to  the  total  power  carried  in 
the  device  for  various  parameters,  is  plotted  as  a  function 
of  frequency  in  Figures  9  and  10. 

Several  facts  can  be  studied  from  the  graphs. 

(a)  The  amount  of  power  carried  in  each  region  is 
frequency  dependent. 

(b)  As  frequency  increases  relatively  more  power  is 
carried  through  region  2,  the  dielectric  substrate. 

(c)  Relatively  more  power  is  carried  through  region  1 
as  the  permittivity  is  lowered. 

(d)  In  the  low  frequency  region,  more  power  is  carried 
through  region  1  as  the  substrate  thickness  increases. 
At  high  frequencies  the  power  carried  in  region  1 
approaches  the  same  value,  for  all  values  of 
substrate  thickness. 

(e)  As  strip  width  increases  more  power  is  carried 
through  region  2,  the  dielectric  substrate. 
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Fiaure  7.  Average  Power  Distribution  in  Region  1  Vs. 

n  (an)  for  f  =  AGHz,  D  =  0.125",  W  =  0.45", 
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Figure  8.  Average  Power  Distribution  in  Region  2  Vs. 

n  (an)  for  f  =  APnHz,  D  =  0.125",  W  =  0.45", 
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I  IK  PERTURBATION  ANALYSIS  OF  EDGE-GUIDED 
MODE  ISOLATOR  ON  FERRITE  SUBSTRATE 


A.  PERTURBATION  EXPRESSION  FOR  PROPAGATION  CONSTANT 

Perturbati onal  analysis  is  used  to  determine  solutions 
for  perturbed  problems  slightly  changed  from  other  problems 
to  which  solutions  are  known.  The  study  that  was  done  up 
to  this  point  obtained  a  solution  for  the  propagation 
constant  of  an  electro-magnetic  wave  traveling  in  a  struc¬ 
ture  as  shown  in  Figure  3.  The  goal  of  this  study  as  was 
stated  in  the  introduction,  was  to  investigate  the  character¬ 
istics  of  waves  traveling  in  the  same  structure  as  shown  in 
Figure  3  but  built  on  a  ferrite  substrate.  This  goal  could 
be  achieved  by  using  the  perturbation  approach  considering 
the  dielectric  case  as  the  unperturbed  problem,  and  the 
ferrite  case  as  the  perturbed  problem. 

A  pe rturbati ona 1  expression  for  the  propagation 
constant  due  to  small  changes  of  a  material  type  in  a 
guidance  structure  is  given  in  reference  8  as 


Y  '  +  Y* 


joi  //  [(e.[Ax-]-E'-E*)  +  (u.[Ax„]-H' -H*)]  ds 

AS  _ _ _ ^ -  (133) 

//  (E*xH'  +  E'xH*)-a.2  ds 


where  primes  denote  perturbed  quantities. 
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When  external  magnetic  susceptibility  is  used,  and 
assuming  that  the  dielectric  constant  remains  unchanged, 
hence  [^x^]  =  0  then  unperturbed  magnetic  and  electric  fields 
are  used  instead  of  the  perturbed  ones,  and  equation  (133) 
can  be  rewritten  as 


ff  ds 

AS 

//(E*xH  +  ExH*)-a^  ds 


(134) 


In  the  above  equation  the  following  notations  are  used 


y'  =  a'  +  j3'  perturbed  propagation  constant 


Y  =  a  +  j  3 
AS 

S 

e 


unperturbed  propagation  constant 

cross  section  of  waveguide  which  is 
perturbed  by  a  change  of  the 
material  type 

waveguide  cross  section 

external  magnetic  susceptibility  tensor 


H  and  E 


unperturbed  magnetic  and  electric 
f i  el ds 


The  external  D.C.  magnetic  field  is  applied  perpendicu¬ 
lar  to  the  direction  of  propagation  as  shown  in  Figure  11. 
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Not«:  tht  upp«r  wall 
of  wavaguidi  wot 
romovfd* 


Figure  11.  Shielded  Isolator  Structure. 
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According  to  reference  1  the  lines  of  the  magnetic 
field  lie  in  Z-X  plane  hence  an  interaction  between  the  R.F. 
field  and  the  ferrj'te  results  when  the  D.C.  magnetic  field 
is  applied  perpendicular  to  Z-X  plane. 

The  direction  of  the  applied  magnetic  field  implies  no 
R.F.  demagnetization  in  the  y  direction.  Thus  the  external 
magnetic  susceptibility  tensor  is  given  by 


[x®]  = 


X  0 

^xx 


'xz 


^zz 


^xxaxax  ^xzaxaz 

+  \f  -  -  +  Y  -  - 

^zxazax  ^zzazaz 


(135) 


and  the  integrand  in  the  numerator  of  equation  (133)  can  be 
explicitly  written  as 

Cx®]-H-H  =  +  XxzH^]H*  +  +  X^z^z^Hz*  (136) 


The  denominator  of  equation  is  simply  4Pave  so  one  can 
rewrite  equation  (134)  as 


Jwy. 


D  A 

0  0^^^ 


(137) 


+  X^zl^zl 


"  x.zHz^:  ^XzxHxH*  )dxdy 


4  Pave 
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B.  MAGNETIC  SUSCEPTIBILITY  TENSOR  ELEMENTS 

According  to  small  signal  approximation  of  the  equation 
of  motion  in  reference  8  two  frequencies  can  be  defined  as 


“x  =  (“o  -  ^  ^  =  %  -  %• 


(138) 

(139) 


where 

iOq  =  yHoy 

=  Y47rM5 

Y  =  1.4*g[MHz/oe]  =0.879*g[Mrad/oe] 
g  =  Lande'  factor  (-  2) 

Hoy  =  applied  D.C.  magnetic  field  in  y  direction  [oe] 
AttM-  =  saturation  magnetization  [Gauss] 

N^,Ny,Nz  =  Demagnetization  factors  (in  this  case 

Nx  =  Nz  =  0,  Ny  =  1) 


Define  the  ellipticity  of  the  normal  modes  of  the 
uniform  precession  [Ref.  8]  as: 


(140) 


and  Kittel  resonance  relation  for  an  ellipsoid  [Refs.  8  and  9] 
as 


“r  '  V“x‘“2 
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The  Landau-Li fshi tz  damping  factor  is  defined  in 
reference  9  as 


a 


(141) 


where  aH  is  the  line  width  of  the  ferrite. 

In  this  study  the  ellipticity  e^  is  equal  to  1  and 
“r  "  ‘^o  "  taking 

X  =  x‘  -  jx“  (142) 

and  knowing  that  [Ref.  8] 


^xz  “  "  ^zx 


(143) 


then  the  susceptibility  elements  can  be  expressed  as 


^  XX  "  ^  zz 


,22.,  22 
A 


^  XX  "  ^  zz 


+  cji)^(l+a^)] 


(A)  •  -  (jt)^(l+a^)] 


I  -  .  _  4  ni 

X  xz  ■  ■  X  zx  -  J - 


2(A) 

m  r 


^xz”"^zx”'^  A 


and 


,  r  2  2,,^  2>t2  ,  .  222 

A  =  L<J|.  *i<)  (i+ct  ;J  +  4(1)^  (D  a 


(144) 

(145) 

(146) 

(147) 

(148) 
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Since  x^x  ^  ^zz  ^xz  ^  "  ^zx  simplify 

the  perturbation  expression  given  in  equation  (137)  as 
f ol 1 ows 

I  ?  Cx,,(|Hxi2  .  |H^|2) 

y'  +  Y*  =  jcoy^  - -  (149) 


C.  COMPUTATION  OF  PERTURBATION  EXPRESSION  IN  THE 
SPECTRAL  DOMAIN 

In  order  to  solve  the  perturbation  expression  [Eq.  149], 
the  magnetic  fields,  and  were  substituted  in  terms  of 
the  scalar  potentials  as  follows 


★  jwy  DAi  9<J>5(x,y) 

"  =  4Pl7e  X,,[{ljS- 

00  ' 


9X 


+  jwe2 


9(|>|(x,y ) 
9y 


)  (150) 


(+j6 


34)2(^»y)  9(|>|(A,y;^  .  ^  4  h/ 


3X 


-  ja)e2  - )  +  KC2  ({>2(x,y)]  + 


+  2x^^[Kc2^*5(x,y)(?jS 


3<i'2(x,y) 

3X 


-  jwe. 


3<('|(x  ,y). 

W 


)]ydxdy.. 


The  upper  and  lower  signs  denote  waves  traveling  in  negative 
and  positive  z-directions  respectively. 

Apply  Parseval's  theorem 


+  Y*  = 


4PaveA 


ri  =  -00  *0  ^ 


(151 ) 


+  j6ana)e2('^2^°^'^*^^' 


9^9(°^n,y)  .*  9$;(an,y) 
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(151  ) 


2  0 
111  £2^ 


3$!(an,y)  ^ 

^35 — I  : 


+ 


J.h, 


2XxzKc|$2(“'’ ,y  )[  +  6an$2*(an  ,y)  -  ju£2-g|-3  }  dy . 


3$; 


3y 


The  integration  with  respect  to  y  was  computed 
analytically  for  both  regions  where  is  either  a  real  or 
an  imaginary  quantity. 

Recall  equations  (49)  through  (52) 


.  (CM(an)  cosh  Ypy 

^2(“n,y)  2 

(c!f (an)  cosY2"y 


Y2  real 
Y2  imaginary 


3 $2  («ri  ,y ) 

ay 


Y2D^(cin)  coshY2y 
Y2"Dy(an)  cosY2"y 


Y2  >^eal 
Y2  imaginary. 


Substituting  into  equation  (151)  for  both  hyperbolic  and 
trigonometric  cases,  and  integrating  with  respect  to  y,  one 
can  obtain: 

For  the  hyperbolic  case 

'  sinh2Y2D 


(152) 


Jwy  \  :>irin^Yoi^  9  9  9  h  ^ 

Z  2y,  "  "  Kc2)|C„(an)|  + 


0  0  9a  ^  h 

j  £2  Y2  I  ^  I  ^ 


-  c5*(an)  D®(an)]  - 


-  j2Xx2^‘^2  “  '^2^2'^H  C^(an)[ 


^  sinh2Y2D 


2y, 


+  D] 
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And  similarly  for  the  trigonometric  case 


(153) 


•  * 
Y  +  Y 


J“^0 

SPaveA 


M 

sin2Y2  D 

-^77- 


+  D][(e^an^ 


+  KC2^). |c{(an)| 


2 


+  2^Y2  I  ^7  ( ^  6ana)e272  (^7(“'^)^7  (“fi) 

+  c!f*(an)D7an))]  - 

2  "  e*  h  Sin2Y2"D  ) 

-  “£2^2  Dj  (cin)C-j-(an)[ — -  +  D]  |  . 

D.  COMPUTATION  OF  NORMALIZED  PHASE  CONSTANT  AND  ATTENUATION 
The  perturbed  and  unperturbed  propagation  constants  have 
the  form  of 

y'  =  ct'  +  jS'  (154) 

Y  =  j3  (a  =  0  for  dielectric  case) 

Thus  the  left  hand  side  of  the  perturbation  expression  is 
equal  to 

y'  +y*=  a'  +j3'  +  (-ja)  =  a'  +j(3'  -  3)  (155) 

Using  equation  (155)  one  can  obtain 

•  ® 

n  =  -co 
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(a) 


FREQ  (GHz) 


Figure  12.  a)  Attenuation  and 

b)  Normalized  Prooagation  Constant  for  Waves 
Travelinq  in  +2  and  -Z  Directions  for: 
ah  =  75  0£„  ^  =  1916.3  oe , 

4nMS  =  1200  Ga,  g  =  1 .99  =  15.2 
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Q  (db/cm) 


Figure  13.  a)  Attenuation  and 

b)  Normalized  Propagation  Constant  for  Waves 
Traveling  in  +Z  and  -Z  Directions  for: 
aH  =  75oe,  Hpi  p  =  2271. 42oe, 

4nMS  =  1200  =  1.99,  =  15.2 
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Figure  14.  a)  Attenuation  and 

b)  Normalized  Propagation  Constant  for  Waves 
Traveling  in  +Z  and  -Z  Directions  for: 
aH  =  75  oe ,  Q  “  2628.6  oe , 

4tMS  =  1  200  Ga,  g  =  1.99,  =  15.2 
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or 


00 


SPave  •  A 


(157) 


n  =  -00 


oo 


(158) 


Detailed  development  of  the  last  two  equations  is  given  in 
appendix  B. 

Figures  12,  13  and  14  were  plotted  for  three  different 
values  of  external  magnetic  fields  Hoy,  hence  for  three 
different  resonant  frequencies. 

Several  facts  can  be  studied  from  these  plots. 

a.  The  ratio  of  backward  to  forward  attennuation  is 
very  large,  as  required  for  an  isolator. 

b.  Both  backward  and  forward  attenuations  increase 
as  the  resonant  frequency  increases. 

c.  The  ratio  of  backward  to  forward  attenuation 
increases  as  the  resonant  frequency  increases. 

d.  The  bandwidth  of  the  isolator  is  determined  by 
the  ferrite  linewidth. 

e.  The  normalized  phase  constants  in  both  directions 
are  frequency  dependent  and  exceed  high  values  as 
approaching  the  lower  side  of  the  isolator's 
bandwi dth  . 
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The  region  of  negative  values  of  normalized  phase 
constants  was  not  investigated,  thus  at  present  no  explana¬ 
tion  is  obtained.  .Hopefully  laboratory  measurements  will 
indicate  whether  it  was  an  error  due  to  the  inherent 
approximations  of  the  perturbati onal  technique,  or  how  large 
are  the  deviations  from  the  correct  results. 
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IV.  COMPUTER  PROGRAMMING 


The  computer  program  that  was  developed  for  this  study 
was  written  in  FORTRAN  IV  language  and  has  three  major  steps: 

1.  Computing  wavelength  ratio  X'/X,  for  waves  traveling 
in  a  device  built  on  dielectric  substrate  (Fig.  3). 

2.  Computing  power  flow  in  a  device  built  on 
dielectric  substrate. 

3.  Computing  normalized  phase  constant  3F/3,  and 
attenuation  in  db/cm  for  waves  traveling  in  a 
device  built  on  ferrite  substrate  (Fig.  3). 

These  computations  were  done  for  both  forward 
and  backward  directions  of  propagation. 

The  computations  of  the  wavelength  ratio  in  step  1  were 
made  by  finding  the  root  3  that  solved  the  de termi nental 
equation  (113).  First,  equation  (113)  was  solved  for  an 
arbitrary  value  of  3,  and  then,  by  the  use  of  Newton-Raphson 
iteration  method,  3  was  computed  until  the  change  between  two 
3's  in  two  successive  iterations  was  in  the  6th  or  7th  digit 
after  the  decimal  point.  From  the  value  of  the  last  3, 

X'/X  could  easily  be  computed. 

In  step  2  the  power  flow  was  computed  for  both  hyperbolic 
and  trigonometric  cases,  depending  whether  was  a  real  or 
imaginary  quantity.  The  computed  3,  and  the  coefficient  ratio 
a2/a^  obtained  in  part  1,  were  used  as  entries  in  part  2. 
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In  step  3  both  power  flow  and  computed  3  were  used  as 
entries.  This  step  was  executed  twice  each  time  for  waves 
traveling  in  both  directions.  The  output  of  this  step  could 
be  obtained  for  only  one  value  of  the  following  parameter  at 
each  run: 


Ferrite  linewidth 
Applied  D.C.  magnetic  field 
Saturation  magnetization  and 
Lande'-g  of  the  ferrite. 


All  three  steps  could  be  calculated  for  different 
values  of  frequencies  strip's  width  and  substrate's  thickness 
in  one  run. 

The  required  input  data  to  the  program  is  given  in 
Appendix  C. 

The  limits  for  all  summations  in  the  program  were  chosen 
by  a  trial  and  error  method,  since  there  were  no  definite 
limits  that  could  be  pointed  out. 

The  upper  and  lower  frequencies  of  operation  are  bounded 
as  follows. 

Lower  frequency  is  bounded  by  the  strip's  width  such  as 


X^max 


W 


(159) 


or 


C 


(160) 
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Approximated  upper  bound  was  applied  according  to 
reference  10.  This  reference  discusses  the  case  of  a 
rectangular  guide  with  dielectric  slab  perpendicular  to  the 
electric  field  as  shown  in  Figure  15. 


‘r, 

H 

J: 

D  < 

-  A  - ► 

^  1 

Figure  15.  Waveguide  Filled  with  Dielectric 
Slab  Perpendicular  to  the 
Electri c  Fie  1 d . 


The  guided  wavelength  is  given  as 


X 


g 


>1  lid  - 

n  er  9 


(161) 


The  cut-off  wavelength  of  this  guide  is  the  one  that  set 
the  denominator  equal  to  zero,  so  from  equation  (161)  one  can 
obtain 


^0  mi  n 


(162) 


and 


max 


C 

^0  min 


(163) 
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V.  CONCLUSIONS 


In  this  study'  a  theoretical  analysis  of  a  model  for  a 
shielded  edge-guided  mode  isolator  was  presented. 

The  analysis  was  based  on  a  complete  solution  for  an 
unperturbed  problem  (the  dielectric  case)  and  then  by  the 
use  of  a  perturbation  technique,  a  solution  was  derived  for 
a  perturbed  problem,  (the  ferrite  case). 

Final  results  indicated  that  above  and  below  the 
resonant  frequency  the  forward  attenuation  is  very  low  - 
negligible,  while  the  reverse  attenuation  is  high.  These 
are  good  regions  in  which  to  operate  as  an  isolator.  One 
result  which  is  still  unexplained  is  that  the  normalized 
phase  constant  is  negative  in  the  frequency  region  above 
the  resonant  frequency. 

A  computer  program  was  developed  for  all  steps  of  the 
theoretical  analysis. 
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APPENDIX  A 


AVERAGE  POWER  FLOW  IN  REGIONS  1  and  2 


Equation  (129)  indicates  the  total  power  flow  in  region 
1.  In  order  to  solve  this  equation  both  A®  and  A*’  should  be 
known  explicitly. 

Recal 1 


A®(  n) 


Kc^^sinhy,  (D-H) 


(A-1) 


ones’  (an)  , 

-  ^ - ^v(an)]  - 7^-4 - rfT-uT  (A-2) 

x'  10^1  Y-j  SI  nhy-]  ( D-H ) 

M^^’^(an,6)d'^(an)  +  M2^ (on  ,  6  ).f2  ^ 

S^(an)  =  M3^’''’(an,B)J-^(an)  +  ’ ''’ ( on  ,  6  )^2  ( on  )  (A-4) 

Since  M^^’^(an,B).  M2^’^(an,B),  M2^’^(an,B)  and 

H  T 

’  (an,6)  are  imaginary  quantities,  one  can  define 


A^an)  = 
^^(an)  = 


(an  , 

II 

( an  . 

II 

CQ 

M  ^ 

*'*'(an  : 

,3)  = 

jin^^’^(an,B) 
jm2^ ’^(ai , B ) 
’^  ( on  ,  B  ) 


(A-5) 


M^^’^(an,B)  =  jm4”’^(an,B) 
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Furthermore,  referring  to  equation  (115)  one  can  assume 
that  a 2  is  a  real  quantity  and  a^  is  imaginary  quantity • 
So  one  can  write  . 


=  ja^  (A-6) 

where  a^  is  real  quantity. 

Substituting  equation  (A-5),  (A-6),  (110)  and  (111) 
into  equations  (A-3)  and  (A-4)  one  can  obtain 


S^(an)  =  5*  (an)  +  jE  (an) 

A  A  -j  A  2 

S’  (an)  =  E  (an)  +  jS  (an) 
z  Z2 


where 


(ctn) 

^1 


5*^  (ctn) 
X2 


(an) 

^1 


a-jC-  —  m^*^*^(an,3)g^Q(an)  - 

"l 

m2^  ’'*'(an  ,3  )  ^  ^ 

•  m^^’^(an,3)  9j^g(otn)  - 

^1 

m2^’^(an  ,3)  f^QCotn)] 

a^[-  m^*^ ’^(an  ,3)  ‘ 

^1 

m^^ ’'^^(an  ,3)  f^gCoin)] 


(A-7) 

(A-8) 


(A-9) 


(A-10) 


(A-11) 
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(A-12) 


s.  (an) 
^2 


,g)g^g(an)  - 
m^^.’^(on,e)f2o(an)] 


and 

w 

y^m/’'''(an.6)B(an) 

^2  _  n=  _ 

i,  A  H  T  2 

’^(an  ,6)  I  g^(an)  I 

n=-«> 


(A-13) 


After  investigating  equations  (A-9)  through  (A-12)  it 

could  be  seen  that  E  (an)  and  E_  (an)  are  even  functions, 

^2 

while  E  (an)  and  E  (an)  were  odd  functions.  Substitution 
•^2 

of  these  functions  into  equations  (A-1)  and  (A-2)  led  to 
explicit  expresions  of  A  (an)  and  A  (an),  hence  the  power 
flow  in  region  1  could  be  calculated. 

It  is  clear  that  the  power  flow  in  region  2  for  both 

the  hyperbolic  and  trigonometric  cases  could  be  calculated 

0  h 

in  the  same  way  since  the  coefficients  D^(an),  Cjj(an), 

o  h 

Dj(an)  and  Cj(an)  are  functions  of  5'^(an)  and  £’^(an)  which 
their  explicit  expressions  are  given  in  equations  (A-9) 
through  (A-12). 
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APPENDIX  B 


NORMALIZED  PHASE  CONSTANT  AND  ATTENUATION 


As  was  stated'in  equation  (157)  and  (158) 


SPaveA 


(B-1) 


00 


^  SPaveA  Z 


(B-2) 


n  =  -00 


So  in  order  to  calculate  a'  and  8'/8>  both  I^  and  should 
be  known  explictly. 

For  the  hyperbolic  case>  equati on  (151)  had  to  be  divided 
into  real  and  imaginary  components. 

Recall  equations  (67)  and  (68)  and  substrate  equations 
(A-7)  and  (A-8)  for  5’^(an)  and  e^{olT])  one  can  obtain 


* _ { r  “  ^  ^ 

cj|i2Y2^^  nhY2^  ^  ^2 


s  (an)  +  E  ^  (ccn)]  +  (B-3) 

Z2  X2 


+  j[-"  (an)  -  E  (an)]} 

KC2  1  1 


[s,  (an)  +  (an)] 


(B-4) 
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Substituting  equations  (B-3)  and  (B-4)  into  equation 
(152)  and  separating  to  real  and  imaginary  parts,  the 
following  equations  are  obtained 

2  2  2 

3  an^  +  Kc/  „  ^  2 

"  -2  2  2  +  (B-5) 

0.  U2  ^2  ^  2  2  "2 


2  2  2 

0)  ep  Yo  9  9 

- 4 -  [^7  (ctn)  +  E  (an)  + 

^1  ^2 


“26an  r3an  pt-j  2/  \  ,  „  2/  \n 

+  - ^  { - p-  [S'  (an)  +  E_  (an)j  - 

KC2  K  2  1  2 


-  S’  (an)S'  (an)  -  Ej  (an)  (an)} 


Zi' 


Rip  =  —  [s’  (an)s  (an)  -  (an)r  (an)] 
^  ^2  ^1  ^2  ^2  ^1 


(B-6) 


ctghypD  .  x'  x' 

RE  =  +  - 5 - ][-|^  RI,  +  ^  Rio]  (B-7) 

^2  sinh^Y2^^  ^  \  d  d 


ctghvpD  n  x" 

Im  =  C - ^  T 

^2  sinh^YoD 


1  ■  “T'  '^^2^ 
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Similarity  for  the  trigonometric  case  the  coefficients 
h  0 

Cj(an)  and  Dj(an)  can  be  expressed  in  terms  of  £'^(an)  and 
5'^(an)  as  follows 


c!l(„n)  =  - ^ -  {[^  E  (an)  -  E  (an)]  +  (B-9) 

a)U2Y2^inYp  Kc2  2  ^2 


+  j[-  ^  E^  (an)  +  (an)]} 


Dj(an ) 


j  2  "  5^_(cin) 

Kc  2  s  i  n  Y  2  D 


(B-10) 


and  both  the  real  and  imaginary  parts  of  equation  (152)  have 
the  form  of 


RI 


1 


.2  „2  ,  .p  2 

S  an  +  KC2 

2  .  2  „2 

w  y2  ^2 


E  (an) 
^2 


2 

(an)]  + 

X2 


(B-11) 


+  [-  ^E^  (an)  +  E^  {an)]S  + 
KC2  1  1 


2  2  "2 

OJ  £  rt  Y  o  9  9 

+  - (an)]  + 

KC2^  ^1  ^2 


-  2_3a^  Ie  ^(an)  +  S'  ^(an)]  - 

Kc2  Kc2  ^1  ^2 

-  S',  (an)S  (an)  -  E  (an)S'  (an)} 

Zi  Z2  X2 
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(B-12) 


RI 


2  y2 


[E  (an)g-^  (an)  -  E^  (a.n)E^  (an)] 


^2 


r  D  ^  ctgyoD  x‘  -Xy-7 

sin^Y2D  ^2  \  d  d 


(B-13) 


SI  n 


n  CtQYoB  Xwy  Xy^ 

^  1  r.  XX  pr  _XZ  p  T  - 

2  ...  ^  2  2  ^4- 


■yP  ^2 


(B-14) 


After  substituting  properly  equations  (B-7),  (B-8),  (B-13) 
and  (B-14)  into  equations  (B-1)  and  (B-2)  the  later  can  be 
solved  for  the  attenuation  and  normalized  phase  constant. 
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APPENDIX  C 


COMPUTER  PROGRAM 


THEORY  OF  FIELD  D ISPL  EMcNT  DEVICES 

SHEILDEO  EDGE  GUIDED  MCCE  ISOLArOR 
(BOX'S  DIMENSION  0.4* **0.9**) 


THIS  PROGRAM  HAS  THREE  MAJOR  STEPS  CORRESPONDING  TO 
THE  STEPS  OF  THE  THEORETICAL  ANALYSIS  OF  THE  DEVICE. 
THE  THREE  STEPS  ARE: 

1.  COMPUTING  WAVELENGTH  RATIO  OF  THE  DEVICE  BUILT  ON 
DIELECTRIC  SUBSTRATE. 

2.  COMPUTING  POWER  FLOW  IN  THE  DEVICE. 

3.  COMPUTING  NORMALIZED  PROPAGATION  CONSTANT  AND 
attenuation  in  CB  per  CM  OF  FERBITE  BUILT  DEVICE 
BY  USING  PERrjReATICN  THEORY  ANC  DATA  CBTAIN'EC  IN 
STEPS  1  G  2. 


PROGRAM  ACCEPTS  FOLLOWING  DATA: 

1. FIRST  CARD-LltMl,MNl  FORM AT ( 3 ( I  2, X2 ) ) 

Li  -NUMBER  OF  FREQ.  CATA  CAROS. 

Ml  -NUMBER  OF  STRIP  WIDTHS  DATA  CARDS. 

MNl-NUMSER  OF  SUBSTRATE  THICKNESS  DATA  CAROS. 

2. SECOND  CARD-(FREQ(L) fL=l,Ll )  F0RMA^(D9.2) 

FREQ(L)  ARE  THE  FREQUENCIES  IN  HERTZ  AT  WHICH  THE 
COMPUTATIONS  ARE  EXECUTED. 

NUMBER  OF  * 'SECOND**  CARDS^Ll. 

3. THIRD  CARO-(W(M) tM=l ,M1)  FORMAT  (DIO. 4 ) 

W(M)  ARE  THE  VARIOUS  STRIP  WIDTHSt  IN 
NUMBER  OF  "THIRD**  CARDS  =  M1. 


4. FOURTH  CARD-(D(MN),MN=1,MN1)  FORMAT (010.4) 

D(MN)  ARE  THE  VARIOUS  SUBSTRATE  THICKNESS,  IN  METERS. 
NUMBER  OF  *» FOURTH**  CARDS^MNl. 


5. FIFTH  CARD-EPSR  ^0 RM AT ( F5 , 2 ) . 

EPSR-  RELATIVE  DIELECTRIC  CONSTANT. 


6. SIXTH  CARO-CELTH  FG RM AT ( D12 . 3  ) 
DELTH-LINEWIDTH  OF  FERRITE  IN  CERSTEDS. 


7.  SEVENTH  CARD-HO  i=ORM  AT  (  D12 . 3  ) 

HO-APPLIED  O.C  MAGNETIC  FIELD  IN  OERSTEDS. 


8. eighth  CARD-AMAGS  F0RMAT(D12.2) 
AMAGS-SATURATION  MAGNETIZATION  IN 


GAUSS • 


9. NINTH  CARD-G  FORMAT (F4. 2) 
G-  LANOE-G  OF  FERRITE. 


OUTPUT : 

1.  WAVELENGTH  RATIO  FOP  DIELECTRIC  SUBSTRATE. 

2.  NORMALIZED  PROPAGATION  CONSTANT  FOR  BACKWARD  AND 
FORWARD  DIRECTIONS  OF  PROPAGATION. 

3.  ATTENUATION  FOR  BACKWARD  AND  FORWARD  DIRECTIONS  OF 
propagation. 


THE  OUTPUT  IS  PRINTED  FOR  ALL  FREC'JENCIES,  STRIP 
WIDTHS  AND  SUBSTRATE  THICKNESS. 

CAUTTON:  THE  UPPER  AND  LOWER  FREQUENCIES  OF  OPERATION 
ARE  BOUNCED  TO  ELIMINATE  UNOESIREC  MODES  GF 
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propagation  in  the  device. 

ALL  COMPUTATIONS  ARE  PERFORMED  IN  CC'JSLE  PRECISION. 
LANGUAGE-FORTRAN 

PROGRAM  DEVELOPED  3Y  LT.J.G  RAM  SHARON  ISRAELI  NAVY. 
THESIS  ADVISOR-  PROF.  CR.  J.3.KNQRR. 

NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CALIFORNIA  92940 
FEBRUARY  1976 

IMPLICIT  REALMS  (A-H,C-Z) 

DIMENSION  T0T(4),FREQ(1Q),W(8)  ,C(6) 

ZZ=20.D0 

S=0.D0 


FORMAT  (D9.3) 
FORMAT  (DIO. 4) 


FORMAT  (3(12, 2X)) 
FORMAT  (F5.2) 

10  FORMAT  (C12.3) 


7  FORMAT  (F4.2) 

READ  (5,3)  L1,M1,MN1 
READ  (5,1)  (FREQ(L) ,L=l,Ll) 

READ  (5,2)  (^i(M) 

READ  (5,2)  (D(MN)  ,MN=1,MM) 

READ  (5,5)  EPSR 
READ  (5,10)  D6LTH 
READ  (5,10)  HO 
READ  (5,10)  AMAGS 
READ  (5,7)  G 

PI=5. 141592654D0 

AMU1  =  4.D0=!''3I^1.  DD-7 

AMUR=1.D0 

AMU2  =  AMU1’!'AMUR 

EPSl=l.D0/(36.00*PI^1.0+9) 

EPS2=EP$1*EPSR 

STEPl-WAVELENGTH  RATIO  CALCULATIONS, 


H-BCX'3  HEIGHT 


H=1.016D-2 


A-BGX*S  WIDTH 
A=2.286D-2 


FPEGUENCY  LOOP 

DO  4000  L=1,L1 
FREQi=FREQ(L)/1.00+9 
0MEGA=2.C0^PI-FREC(L) 
aM£GAl=QMEGA/1.00+9 
WRITE  (6,4)  , 

4  FORMAT  (//////////) 

PRINT  OUT  OF  INPUT  DATA 


199 

200 
2C1 
202 
2C2 
207 


WRITE  (6,199) 

FHRMAT 

MklTB  (67200)  FREQl 

FORMAT  (20X,aH^  FR  cQl  =  ,  F5 . 3 , 7H  ^3HZ,2H 
WRITE  (6,201)  CMEGAl 

FORMAT  (30X,9H=!‘  QMcGA  1=  ,  F6  •  3 ,  5H  GRPirZH  =« 
WRITE  (6,202)  H 

FORMAT  (30X,4H*  H=,E11,4,5H  MTR,2H  -) 

FG^AT  ^  (  30x74H*  A  =  ,E11.4,5H  MTR,2H  -) 
WRITE  (6,207)  EPSR 

FORMAT  (20X,7H=!'  EPSR=,F4.1,9X,2H  ^) 


^  ) 
) 
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WRITE  (6,206)  DELTH 

2C6  FORMAT  (30X,i0H=}'  0ELTH=  ,F7.2,5H  CE  *) 

WRITE  (6,210  HO 

210  format  (30X,10H=!‘  C.C  •  FL0=  ,  F7. 2 , 5H  OE 
WRITE  (6,212)  AMAGS 

212  FORMAT  (30X,10H’i«  S  AT  .  M  AG  = ,  F7 . 2 , 5h  GA 
WRITE  (6,213)  G 

213  format  (3CX,10H*  LAN0E-G  =  ,F4. 2 , 7X ,1H*  ) 
weiTE  (6,208) 

208  FORMAT  (  30X  ,  ) 

SET  HIGHER  FREQUENCY  OF  OPERATION 

GM  =  2.D0*A/DSQRT(  1 .  DO-0  ( MN ) /H’!' (  l.DC-1.  DC/EPSR)  ) 
FREQU=3.D8/GM 

IF  (FREQ(L).GE.FPEQU)  GC  TO  400C 
STRIP  width  loop 
DO  5000  M=1,M1 

SET  LOWER  FREQUENCY  CF  OPERATION 

FREQL=(3.0C+08/DSCnT(5PSP) )/(4^W(M)  ) 

IF  (FREQ(L).LE. FREQL)  GC  TO  5000 
WRITE  (6,205)  W(M) 

205  FORMAT  ( lOX , 2HW  =  , FI 1 . 5/ / ) 

SUBSTRATE  THICKNESS  LOOP 


DO  6C00  MN=ltMNl 
WRITE  (6,209) 

2  09  FORMAT  (  20X  ,  14HL  AMBD  A  *  /  LAMBDA  ,  "'X  ,  IHC  ,  15 X  ,  4hBETA/  ) 


Nl=101 

11=201 


beta-prop AGA7I0N  CONSTANT* 


C 

lOCO 

C 


600 

700 


C 


16 

9 


C 


6 


QBETA=1.0D-1 

eETAl=300.D0 

TOT( 1 )=20.D0 
DO  2000  K=l,2 


BETA=3ETA1 

IF  (K.EQ.2)  GO  TO  600 
GO  TO  700 
eETA=3ETA+08ETA 
AK1=0MEGA=^=0SGRT(  AMU1’»=EPS1) 
AK2=QMEGA^0SQRT( AMU2^EPS2) 
AKC13  =  AK1^^2-6ETA=«'’f'2 
AKC2S  =  AK2*^2-6cTAs;^>;i2 

SUM1=0«000 

SUM2=0.000 

SUM3=0.D00 

TT?WR1=0.D0 

TTPWR2=0.00 

TGTPWR=0.D0 

GO  TC  9 


Nl  =  51 
11=101 
RF11=0.00 
RE22=0.00 
ZM11=0.00 
ZM21=0.D0 
.N  =  -N1 

DO  100  1=1,11 
N=M+1 

ALFAN  =  2  .DO^^PI^N/ A 

GAMA  1  =  DSQRT  (  ALFAN=<^^2-  AKCIS  ) 
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GAMA25  =  ALFAN*={=2-AKC2S 
GAMA2  =  0SQRT (DABS ( GAMA2S  )  ) 

A1  =  1.D0/DTANH(GA.MA1^(H-C(MN)  )  ) 

A2=0‘1£GA*AMUl=i‘GAMAl 

A3=OMEGA^A:'iU2*GAMA2 

A4=ALFAN*BcTA 

A5=0MEGA=^'E?S1^GAMA1/AKC1S 

A6=O.M£GA^£?S2^GAMA2/AKC2S 

IF  (GAMA2S«GE.O,OCCJ  GO  TO  20 

TRIGONOMETRIC  CASE 

A7=DTAN(GAMA2=?'D(Mi\)  ) 

F1=-AKC  1S^A1/A2  +  AKC2S/  (  AT’i'AB) 
F2=A4=!=Al/A2-A4/(A7t^A3  ) 

F3=-F2 

F4  =  A4’S'^2^Ai  /  (  A2'^AKC1S  ) -A5’}'Ai-A4=J'*2/ (  AKC2S^A7*A3  )■ 
GO  TO  30 

hyperbolic  case 

20  A8=1.00/DTANH(GAMA2+0(MN) ) 

F1  =  -AKC1S^A1/A2-AKC2S=^A8/A3 
F2  =  A4^A1/A2  +  A4=^''A8/A3 
F3=-F2 

F4=A4*=!'2*A1/  (  A2^AKC1S  ) -AS^^' A1  + A4=J'=i' 2^A3/ (  AKC2S^*A  3  ) - 
30  DN=F1*F4-F2’pF3 
AM1=F4/0N 
AM2=-F2/0N 
AM3=-AM2 
AM4=F1/DN 

SN  =  OSIN(ALFAN’<^W(M)  ) 

CN=DC  JS  (ALFAN^i'W  (M  )  ) 

R1=DEXP(ZZ) 

R2=(  ZZ/W(M)  J^i‘2  +  ALFAN’»'^2 

R3=R1=^CN 

R7=R1*5N 

Q=(  (PI/(2.D0=^'W(M)  )  )=^‘=^'2-ALFAN’!'=^'2) 

GX1  =  PI/  (2.D0=«'W(M)  j^i'CN/G 
GX2=( PI/(2.00^W(M) )^SN+ALFAN) /G 
FZ1=(  ZZ/W(M)  =^‘(R3-1.00)  +ALFAN*R7)  /R2 
FZ2  =  (  ~ALFAN=i=(  R3-1  .DO  ) +ZZ/W(  M  )  =i'R7  )  /R2 
IF  (0ABS(T3T(1)  I  .G6.1,D+01)  GO  tq  a 

STEP  2  -POWER  FLOW  CALCULATIONS. 

AMP1=1.D-10 
Z1=-SUM2/SJM3 
ZS^Zl’^GXl 
Z2  =  Z1=!‘GX2 

EZ1=AMP1^(  -AM3’!'Z2-AM4^FZ1  ) 
EZ2=AMP1-(A^3=^'Z3-AM4^FZ2) 

EX1=AMPL*( -Z2^AM1-AM2*FZ1 ) 

EX2=AMP1^(Z3«AM1-AM2*  FZ2) 

EZS  =  EZl*=i'2+£Z2’5‘*2 
eXS=EXl=^=i‘2+EX2--2 
IF  (S.EQ.l.DO)  GO  TO  13 

REGION  1  FREE  SPACE 

V0=GAMA1*(H-C(MN) ) 

V1  =  0MEGA=*6ETA^EPS1/2.D0 
V2=1.DO/(OTANH(  VO  >=i'GAMAl) 

IP  (VO. GE. 85. CO)  GO  TO  12 
V3  =  (H-D(m.N)  )/((DSINH(V0)  )^^2) 

GO  TO  11 
12  V3=0.Da 

11  V5  =  B6TA/(2.D0-0MEGA=i«Af^Ul) 

V4=ALF  AN=«'eETA/AKC  IS 

V6=(  V4=i‘EZ2-EX2)  -*2+(-VA^EZl  +  EXl  )-  +  2 
V7  =  V4^EZS-EZ1=fEX1-EZ2^EX2 
TR11=-V1=«'ALFAN*=^'2/(AKC1S''J'^2)=^‘EZS*  ( V2-V3  ) 

TR12  =  -V5=i'ALFAN^-2/(GAMAl=<'^2)-V6’l‘(  V2+V3) 
^Ri3=-V5’«'V6=F(  V2-V3) 


A6/A7 


A6=i'A8 
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14=-V1-(GAMA1=^*=!'2)  /  (  UKCIS^'^2  )  -F Z S*  ( V2+ V3  ) 

156=  ALFAN/(0M£GA’«'AMUl=i‘AKC15  )  ’i' (  3  ETA'^=’!^  2  +  A  K  2  )  ^VT^S'VZ 

R1  =  -(TR1  1+TR12  +  TR13  +  T;R  14  +  75  156  )  /  (2*00=S'A  ) 

in.\HK'3C  1C  r\  /'n  to 


TR 

TR.. 

PWRl  .  _ _ _ 

IF  (GAMA2S.LE.0 


.00)  GO  TC  300 


REGION  2  HYPERBOLIC  CASE 

V10=GA.MA2*0(^N) 

Vll  =  0MFGA*BETA=<'EPS2/2  .CO 
V12=l.  00/(0 TANH  (V 10  )  ^^'GA^-AZ) 

IF  (VIO.GE.IJO. DO  )  GO  TC  14 
V13  =  0(NN)/( (DSINH(Vi0  )  J^^2) 

GO  TO  15 

14  V13=0.D0 

15  V15=BETA/(2.00^GMeGA*A^^U2) 

V14=ALFAN*6eTA/AKC2S 

V16=(  V14^'EZ2~EX2)=f'*2+(-Vl4=i'EZ  1+EX  1)=^'=S‘2 
V17=(  Vi4=!=EZS-EZl^£Xl-EZ2^EX2) 
TR21=-V11=<‘ALFAN*^2*EZS/(AKC2S^^2)  ^(  V12-V13) 
TR22=-V15^ALFAN’^‘*2/GAMA2S^Vi6''!'(  V12  +  V13  ) 

TR23=-V15^V16^( VI 2- VI 3) 

TR24=-V11*GAMA2S/  ( AKC  2S«=5'2  )  « (  V  12  + V  1 3 )  Z S 

TR256  =  ALFAN/(0.'^1EGA’!^AMU2’!‘AKC2S  )  ^  (  3  ETA’i'*2  +  AK2=^=J'2  )  ’S'V  IT^J'V  12 

GO  TO  801 

REGION  2  trigonometric  CASE 

3  00  V25  =  l.D0/(  AKC2S=J‘0SIN(GAMA2*D(  MN  )  )  ) 
V20=CMEGA=«'BETA^E!^S2/2.DC 
V21=GMEGA’^BETA^AMU2/2.D0 
DE71  =  V25=J'EZ2 
CET2=-V25^EZ1 
0ETS  =  0ETl’!'-2  +  DET2=S'=^‘2 

V26=ALFAN=s^BETA=!^V2  5/  {0.MEGA’!=GAMA2*A^U2  ) 
V27=1.00/(OMEGA=^AMU2^GAMA2*0S  IN(GAMA2^0(.MN)  )  ) 

CHTl  =  -V2  6>i'EZ2+V27=i'EX2 
CHT2=V26=J'EZ1-V27^EX1 
CHT$  =  CHTl*=i‘2+CHT2^''J*2 

V23=OSIN(2.DO^GAMA2^D(MN)  )/(2.DC*GA^^A2  ) 

TR21  =  -ALFAN*’!''2^V20*0ET5=«'(D(MN  )  -V2a  ) 

TR22  =  -ALFAN*=«'2=*'V21^CHTS1'(0(MN  )  +  V2a  ) 

TR25=-  (GA^^A2=^'^2)=!'V21=!^CHTS’^'(D(  MN  )-V28) 
tR24=-(GAMA2*’5'2)  =?'V20’!'CETS-(0(  MN)+V2e) 

TR25t>=  (B£TA=?'=i=2+AK2-*2)*ALFAN/2.C0=>(0ETl*CHTl  +  0ET2-CHT2  ) 
SCI  PV^R2  =  -(TR21  +  TR22  +  TR2  3+TR24+TR2  56)  /  (2.00-A  ) 

TCTAL  POWER  FLOW  IN  FREE  SPACE  REGION 


C 

C 


C 


TTPWP1=TTPWR1+PWR1 

TCTAL  POWER  FLOW  IN  CIELECTRIC  PEGICN 
TTPWR2=TTPWR2+PWR2 
TOTAL  POWER  FLOW 

TOTPWR  =  TCTPV\R  +  PWP  1  +  PWP2 
GO  TO  100 


8  FZS=(51**2-2.D0*R3+1.C0)/R2 

GXS=(  (PI/(2.00=<‘W(M)  )  )*’!‘2  +  ALFAN=«^’{‘2+ALFAM^DI/Vv(M):i^CS:N 
X(  ALFAN=i'W(M)  )  )/Q^’;'2 
T1N=AM1^GXS 
T2N=AM4*FZS 
6=-FZ I^GXZ+FIZ^GXI 
T5N=AM2’J'B 
SUMi=SUMl+T IN 
SUM2=SUM2+T2N 
SUM2=SUM3+T3N 
100  CONTINUE 


IF  (S. 
IF  (TO 


EQ.l.CO)  GO  -^0  19 
TPWR.NE.O.DO )  GC  TC  17 
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SUM3SQ=  (SUM3)  =*'=^'2 
TOT(K  )  =SUM1=4:SUM2+SUM3SC 
It  iQABS(T0T( 1) ) .LE.l.C+Cl)  GC  TO  3000 
.2000  CGNTINUE 

DTGT=( rOT(2)-TGT( 1) ) /DEETA 
eETAl  =  DAeS(BETAl-TGT (  1)/DT3T) 

60  IP  (BETAl.LE.AKl  )  GO  TC  50 
GO  TO  70 

50  BETA1=BETA1+1.D0 
GO  TO  60 
70  GO  TC  1000 
30C0  V^VLNG=3.0D+8/FREg(L  ) 

RATI0=2.00^PI /(SETAl^WVLNG) 

V^RITE  (6,211)  RATIO, D(MN), BETA 
211  FORMAT  (24X,F6.4,7X,E11.5,7X,E11.5> 

GO  TO  16 
17  S=1.D0 

SIGN=-1.D0  , PROPAGATION  IN  EACKWARC  DIRECTION  (+Z) 

SIGN=-1.D0 
8CC0  Ni=lCl 
11=201 
GC  TO  9 
21  S=0.D0 
6000  CONTINUE 
50C0  CONTINUE 
4000  CONTINUE 


GO  TO  8001 

STEP3-  CALCULATIONS  OF  NORMALIZED  PRGFAGATICN  CONSTANTS 
ANC  ATTENUATIONS. 


18  V14=BETA^ALFAN/AKC2S 

V16=(  V14*EZ2-eX2)  =«‘^2+  (-VlA^i'EZl  +  EXl) 
V17=(V14*EZS-E21=^EX1-EZ2=^EX2) 

V30  =  (  BETA* AL FAN  )=^'^2  +  AKC 25=^*  =5*2 
V31=(0MEGA*AMU2-GAMA2) 
V22=(OMEGA*EPS2=!*G  AMA2  )  ^^2 
V33=V30/V31=*'V16+V32/AKC2S=5'*2'^E2S 
V34=:(-SIGN)  *2.D0=«'V14^EPS2/AMU2=!^V17 
V45  =  EPS2/AMU2^( EX  1* E Z 2- EX2-EZ 1 ) 

IF  (GAMA2S. LE.0.00)  GC  TO  22 


HYPERBOLIC  CASE 


V10=GAMA2^0(^^N) 

V12=1.D0/(DTANH(V10 )*GAMA2) 

IF  ( V10.GE.170.D0  )  GC  TC  24 
V13  =  D(MN)/(  (DS!NH(V10)  )=<-=<=2) 

GO  TC  25 

24  V13=0.00 

25  RE1=(V33+V34)*( V12+V13 J/2.D0 
ZM1=-RE1 

RE3=V45^( V12+V13) 

ZM2=-RE3 
GO  TC  23 

TRIGONOMETRIC  CASE 

2  2  V40=D(MN)/(0S!N(GAma2-C(MN) 

V41  =  1.D0/(DTA.N(GAMA2=«'C(MN)  )*GAMA2  ) 
RE1=(  V33-V34)=<*(  V4C+V41  )/2.D0 
ZM1=-RE1 

RE3=-V45* ( V40+V41 ) 

ZM2=-RE3 

23  RE11=RE11+RE1 
RE22=RE22+RE3 
ZM11=ZM11+ZM1 
ZM21=ZM21+ZM2 
GO  TC  100 
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19  G1=2.3D6*PI^G 
CMGAM  =  Gl=i'AMAGS 
GMGA0=G1*H0 
OMGAR=OMGAO-CMGAM 
CAyiP=Gl’!‘DeLTH/(  2.00=4^0 EGA) 

DELT=(  0MGAR^«2-ayeGA=«^*2*(  1.D0  +  DAMF’«'«2  )  )^^Z  + 
X4.D0*(0MGAR*0MEGA=!'CAyP  )  ^^2 

C  CALCULATI3M  OF  SUSCEPTIBILITY  TEMSQF'S  COMPONENTS 

XXX1  =  0MGAM*0MGAR=<=  (  (OMGAR=!'^2-CMFGA=«'^2  ) +  (  CMEGA’J'OA  WP  )  ^^2  ) 
XXX2=0MGAM*QMEGA=<'-CAMP*  (CMGAR«*2  +  QyEGA^^2=^‘  (1.0C  +  CAMP^=^'2)  ) 
XXZl=GM£GA'i‘OMGAM^(OMGAR=^‘’^'2-OMEGA^:^'2*(  1 .  DO  +  DAM  P*=i=  2  )  ) 
XXZ2-2.DO=i'OMEGA^^2*aMGAy=5'OAMP=5'CMGAR 
RE  =  XXX1^RE11  +  RE22’<‘XXZ1 
ZM=XXX2=<‘ZMll  +  ZM2i*XXZ2 
904  FORMAT  (///) 

WRITE  (6,904) 

WRITE  (6,904) 

V50=0MEGA«AMUl/(4.C0=^'T0TPWR=i'0ELT«A) 

ATTEN  =  -ZM=<'V50 

DB=8.685889638D-02*ATTEN 

BETAP=BETA+RE-V50 

RAT=BETAP/8ETA 

IF  (SIGN.EQ.l.DO)  GO  TO  81 

WRITE  (6,898)  DB 

393  FORMAT  ( lOX , 2 1H3ACKW AP C  ATTENU AT  I CN  =  , El  1 . 5 , 6H  CB/CM//) 
WRITE  (6,897)  PAT 

397  FORMAT  ( lOX , 21H3ET AF/ BETA  3ACKWAPC  =,510.4//) 


C  SIGN=+1.00  , PROPAGATION  IN  FCRWARO  CIPECTION 
C 


(-Z) 


81 

900 

899 

3001 


SIGN=+1.D0 
GO  TO  8000 
WRITE  (6,900)  DB 

FORMAT  (  10X,20HFQRWAR0  ATTEN'JATIGN  = 
WRITE  (6,399)  RAT 

FORMAT  (  10X,20HBETAF/BETA  FORWAPC  = 
GO  TO  21 


STOP 

END 


511.5, 6H  DB/CM//) 
E10.4//) 
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